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Ie INTRODUCTION 


The maximum power level of a nuclear reactor is fre- 
quently limited by the rate of heat removal. The actual de- 
sign of a power reactor represents a compromise between tne 
requirements of heat transfer and those of nuclear physics. 
Wuclear reactor theory is not sufficiently accurate to jguare- 
antee that a particular reactor (especially one of a new 
design) will go critical. Therefore, the nuclear characteris- 
tics of a reactor design must be obtained before construction 
of the full seale reactor. Some of this information can be 
obtained from the use of a subcritical assembly. 

Experiments with uranium graphite suberitical assemblies 
can be carried out to determine the various nuclear constants 
of a proposed reactor lattice. These lattice constanta are 
the material buckling, multiplication constant, lattice dif- 
fusion length, and thermal utilization of the unit cell in- 
cluding the process tube assembly. 

In this investigation six lattice configurations were 
considered, and the buckling of each was determined experi- 
mentally. lattice cells of 6 ine, &5 ine, and l2 in., de- 
noted I, II, and IIil respectively, constituted the three ceo- 
metric lattice arrangements. For each of these, both the 
"wet" and “dry" cases were considered. The measurements and 


calculations with water present in the cooling annulus were 


4 
<a os 













Wee a Few et ee meme qe 
OF TY mm See el eee owed: 9 dba 
——— ee © Oe oo ty Gee) Gems © to ' 
ee i oe : 
i ee 

We © me eethamee| mmdemes qolmapeny © tangy 

a, een oe some Jfeekiios @ (us tied 
STO ltee nS le at Ame elem erm « te onl 
ert Ce. edbwrT Afews ftw) od 
ee ee er 
Nee Hee BONS es ee ee Oe Fem pe mee 
Gi ee tes gutter os! ete as tee ele 
Se eee MEE meiner) Sa ey fe 
Se NOE) PY LE Lhe oe . 
ee ee a 
tt. eT om ON 

WIE else dl ee Chie fee eee Gee: 
Aes ere ue ae 8 84 ed 
— ,20¢5 ‘me 1.8 ww tS ot hw ema Othe 7 

Ae Oe  eeetteee ,prmrmrenrt 10 ame ott ot 

hm li ~ eet -:eeern areal © 
en Ft Owe Oe reas “YS Bee “ee 
Ver lee atte Bf am wee O6ds emehbelete~ 


























a. 


2 


designated as "wet", and those with the process tube, but 
without water, as "dry". 

These six lattice configurations were also considered 
in the theoretical analysis. The buckling of eaeh of the 
six lattice configurations was determined theoretically by 
two methods. ‘“mly one size fuel element and one size proc- 
ess tubing was used in the experimental and theoretical 
analyses of this investigation. 

The methods of Murray (1) and of Rumsey and Volkoff (2), 
(3) were used to determine theoretically the thermal utilize- 
tion and the resonance escape probability. Murray's method 
was shown to give less conservative results than the method 
of Fumsey and Volkoff. The theoretical multiplication con- 
atant derived for each configuration considered was based 
on literature and calculated values of the various diffusion 
lengths, the Yermi ace, and the various macroscopic cress 
sections. The product of ” and © was assumed to be es-= 
sentially constant for the various configurations investi- 
gated. 

Analyses of the data for the given configurations were 


carried out by use of the classical formula for the buckling 


B- = 6 = (1) 
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A comparison between the “wet” and "dry" eases of the 
three lattice cells considered was presented using the experi- 
mental and theoretical values of the buckling. By an analysis 
of the results obtained in this investigation, the effect of 
the water coolent on the nuclear properties of the various 


lattice configurations was determined. 
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II. REVIEW OF LITERATURE 


The theoretical calculation of the multiplication factor 
for a heterogencous reactor with coolant is an extension to 
the method of calculating the thermal utilization and reso- 
nanee escape probability factors for an infinite heterogeneous 
uranium-graphite reactor presented by Weinberg (kh), (5). 
Guggenheim and Pryce (6) extended the method to include a non- 
absorbing coolant and a thin sheath of absorbant aluminum 
around the fuel element. Houston (7) extended the method to 
include neutron absorption in a coolant. Houston's results 
are more suitable for gas or heevy liquid metal cooled re- 
actors due to his assumptions concerning the coolant's negli- 
gible moderating properties and small volume fraction of the 
lattice cell. Rumsey and Volkoff (2), (3) extended the above 
methods further to include the moderating effect of a water 
anmulus. 

Murray (1) presented a simplified method of calculating 
the thermal utilization which considers all but the uranium 
and the graphite in the lattice cell as poisons which reduce 
the basic urenium graphite thermal utilization. Murray's 
method for the calculation of the resonance ascape probability 
does not consider the so called poisons. The Rumsey and Vol- 
koff method appears to be at the present time the most com- 
prehensive theoretical analysis of a water cooled uranium- 


graphite lattice cell. 
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The experiments with subcritical assemblies of the Han- 
ferd Atomic Products Operation included an effort to determine 
the effect of a water annulus around the fuel element on the 
buckling. The effect of aluminum, air, enc other inpurities 
in the lattice upon the buckling was also investigated (8). 

A few measurements were teken using internally socled cylin- 
@rical fuel rods, but the majority of the coolant simulation 
experiments were with solid eylindrical fuel rods of natural 
uranium surrounded by an annulus of water, The objectives of 
the subcritical essembly experiments at Hanford (3) were to 
Getermine experimentally the various lattice constants. The 
existing lattice theory was supplemented and improved where 
possible. Lattice measurements were taken with three dirffer- 
ent fuel rod diameters: 9.925 in., 1.175 ins, and 1.36 in. 
The lattice spacing and the water-aluminum-uranium ratios 
were varied Tor each slug size. 

The results of the series of tests using both "wet" and 
ary" Lattices were canpiled by Clayton (9). A comparison of 
theory and experimental results for one fuel element diameter, 
six lattice spacings, and the "wet* and "dry" configurations 
for the Hanford Atomic Products Operation was done by Gast 
and others (2). Gast used the method of Rumsey end Volkoff to 
calculate the theoretical values of the thermal utilisation 
and the resonance escape probability. 

Related conclusions (8) from the Hanford work with water 


cooled uranium graphite lattices were: 
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(a) For a given natural urenium fuel assembly, 
there 1s a maximum in the buckling vs. C/U 
atom retio curve. The slug dismeter and the 
water annulus thickness are inversely pro- 
portional to the value of C/U at which the 
maximaim occurs. This maximum occurs between 
G/U values of 50 and 100. 

(oh) The buckling "dry" is greater than the buck- 
ling "wet" for larger lattice spacings (c/U 
atom ratio). As the lattice spacing is de~ 
creased, the buckling curves cross and the 
"wet" buckling becomes greater than the "dry" 
buekling. For this lower c/U atom ratio 
region, externally cooled slugs exhibit a 
"fail-safe" behavior on the loss of waber. 

The following conclusions were a result of the Hanford 
experiments(2). The state of the heterogeneous graphite 
uranium lattice theory at present is not completely satis~ 
factory. Hore experimental measurements of individual paran- 
eters such as f and p aro necessary in order to establish the 
basis for an improved theoretical analysis. Corrections ror 
variations in reactor material densities and purities, and 
allowances for holes, such as for control rods, need further 
investigation (2). Finally, the treatment of the energy spec- 
trum of neutrons, fast to thermal, requires more securate 


development of neutron resonance and capture reactions. 
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III. THMORSTICAL AMALYSIS 
Ae General Procedure 


fhe product of the fast fission factor, € , the primary 
fission neutron factor, " , the thermal utilization factor, 
f, and the resonance escape factor, p, is the infinite multi- 


plication constant, ky. 
cyte = 1 6 (2) 


A theoretical analysis of each of the factors would have 
been extremely complicated without certain sisiplifying as- 
sumptions. Thess assumptions are stated where applicable in 
the following theory. Each of the factors was considered 
separately. The product ey was determined and assumec con- 
stant for all the lattice configurations considered. Expres- 
sions for f and p were determined using the methode of Murray, 
Genoted Method A, and the methods of Rumsey and Volkoff, de- 
noted Method B, for tne “wet” and "dry" cases of Lattices I, 
iI and IIi. 

Using values from the literature for the lattice diffu- 
Sion length end the Fermi age with the theoretically deter- 
mined infinite multiplication constant, the buckling was deter- 


mined from Fquation 1, The theoretical buckling was determined 
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for Lattices I, II, and III for both the “wet” and "dry" 
eases using both Methods A and BE. 

The effect of the air annulus and air holes on the neu- 
tron density in all lattice geometries was neglected due to 
the magnitude of the effect in comparison with other material 
effects in the lattice. 

Subscripts used in this analysis were defined as follows: 

us natural uranium 


al slug can, 25 aluminum 


w water 

& graphite 

m moderator 

p process tube, 615 aluninun 
a absorption 

8 scatterings 


re resonance 
sa scattering "absorption" 
f£ slow neutron fission 
t transport theory 
Any deviation of notation from the above list was individuale- 


ly defined. 


B. Fast Fission Factor 


The fast fission factor, € , is the ratio of the total 


number of neutrons produced by fission to the number of neutrons 
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produced by fission by thermal neutrons. Assuming that the 


primary neutron source is distributed approximately uniformly 


over the fuel, a simplified version of the equation given by 


Glasstone 
is 
eass 
where ? 
Te 
Te 
Oy 
and ove 


. + ™ p 
l-pP (i=) 1 oo (no, + o,) 


and Edlund (5, p. 278) for the fast fission faetor 


On 


ge |Z » (y= 2) SP 
C 


c bing Oe JOE 


(3) 
T¢ 


is the probability that ea fission neutron born in 

a rod will make a collision inside the rod in which 
{t was created; 

is the absorption cross section for fast neutrons, 
that is, the sum of the fission and the capture cross 
sections; 

is the elastic cross section for fast neutrons; 

is the inelastic cross section for fast neutrons; 

is the total cross section for fast neutrons. 


o. = Tg, * G6, + DG. 


Using fast neutron cross section values as given by (ug- 


genheim and Fryece (6, p.» 51), Equation 3 reduces to 


0.098 P 
€- b= ae (4) 
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These numerical coefficients are such that when the fuel ele- 
mont radius is less than 1.7 om the caleulated value of € 
agrees ressonably woll with experimental results obtained by 
Argonne National Laboratory (6, pe51). Por natural uranium 
of density 19 gem? and a rod radius of 0.5 in., the value of 
HO, r,, is 0.262, where N is the number of uranium atoms per 
unit volume, and ‘. is the radius of the fuel rod. The quan- 
tity P is a function of TO, . 
uO, r,, shown by Guggenheim and Pryce (6, pe 51) and the above 
computed value of NO, r., the value of P determined was G25. 


- From a relationship of F and 


Solving Foustion for € gave a value of 1.027. This 
value for the fast fission factor was assumed constant for all 


configurations investigated. 
C. Primery Pission Heutrons 
The number of flasion neutrons released per thermal neu- 


tron captured in the fuel is the factor, He Glagastone and 


Edlund define h (5S, pe 83) as 





where “ is the average number of neutrons produced per thermal 


neutron fission; 
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and ae is the total cross section for absorption of thermal 
neutrons. 

The microscopic thermal fission cross section for naturel 
uranium is 4.18 barns. Applying the maxwellian distribution 
correction of 1.128 and the "not - 1/v" factor of 6.961 for 
yas to this value gives an effective fission cross section of 
3.635 barns for natural uranium. The microscopic thermal ab- 
sorption cross section for natural uranium is 7.68 barns. Ap=- 
plying the maxwellien distribution correction and a “not - 1/v" 
factor of 0.990 for natural uranium to this value determines 
an effective absorption cross section of 6.7) barns for natu- 
ral ursnium (1, p. 32). 

Murray (1, pe 95) gives a value of 2.6 for ¥ . Sub- 
stituting the above quantities in Equation 5, the calculated 
value of yy was 1.327. 

It is normal practice in working with subcritical assem- 
blies to experimentally determine k. and then divide k, by 
the product éepf to arrive at a value of " ° The value of a 
found by subcritical assembly experiments was 1.308 neutrons 
per thermal neutron captured (106, p. 85). This value is in 
reasonable arreement to the computed value of 1.327. The val- 
ue of 4 used in the following analysis was 1.368 an@ was as- 


sumed constant for all configurations investigated. 
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D. Thermal UJtilizcetion 


In heterogeneous lattice theory the finite lattice re- 
actor is assumed to be an infinite lattice array. The prob- 
lem of calculating f and p 1s simplified by replacing the 
square colls by equivalent cylindrical umit cells. Table 1 
gives the equivalent cell dimensions for Lattices I, If, and 
ITI. Assuming that little change in the overall neutron 
flux is experienced in traversing one cell, an individual 
cell may be representative of tho total. 

The methods of Murray (1) and of Fumsey and Volkeff (2), 
(3) were used to determine theoretically the thermal utilize- 
tion. ‘These methods were considered seperately and were de- 
noted Methods A and © respectively. 
le Method A 

The depression of the thermal neutron flux in the fuel 
and adjacent moderator complicates the calculation of the 
thermal utilisation. <Averare flux and volume weighting fac- 
tors must both be included in the fractional absorption of 
the fuel. Prom the spatial varietion of the flux in a cell, 
the averare values of the flux cen be determined. 

The thermal neutron diffusion equations for the uranium 


and the moderator from diffusion theory are 


2 
D, Vo By > By 2, + 5, = 0 (6) 
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and 


2 
Py V~ By > Bq 2m + Sq = 9 (7) 


Equations 6 and 7 can be simplified with the following as- 
sumptions: 
(a) The variation of neutron flux at a given 
coll radius willbde the same everywhere. 
(b>) There are no thermal neutrons produced 
in the uranium. Sy is zero. 
(c) The production rate of thermal neutrons 
in the moderator is independent of 
position. Ss is constant. 
(4) The flux is assumed not to vary along 
the fuel cell axis. 
Using the above assumptions and cylindrical coordinates, 


Equation 6 for the fuel becomes 





2 
a“g 1 (% 
—r *s a) ~#,*8 (8) 


where x = f Fl and f is the inverse diffusion length of 


uraniwn, and Foauation 7 for the moderator becomes 





2 - 
- Pn (on - ey — a ) 
ar< r or m 7m Ds (9 





a 











" rege antag 

Magthielpwwemnmmces¥ 
A 
eres me Ue metho tie 
SS Gee ieee oe me ee i 






+ A, a Te ont ot Lh eee oh 
nes om ecemm cay ext © GidAweT lene 
ee ee 


ibe S-+ ay - Gee 


i 








Ly 


where Da ia the diffusion coefficient for the moderator. 
e < ie the square of the inverse diffusion length. 
Equation 8 for the fuel is the modified Beesel equation 
of zero order. FEquation 9 for the moderator is the inhomo~ 
geneous type of the modified Bessel equation of zero orcer. 


The solutions of Fquations @ and 9 are 


P(r) = AI, (# Lr) (10) 
. 2, | 
and pir) = € Ky (a 3) + 6 I, (4 ) *<. (11) 


respectively. The boundary conditions used in the evaluation 
of the constants A, C, and G are neutron flux and current 
continuity et r=r, (uranium rod radius) and the condition 
thet dg _/er = Cat r = rs (the equivalent outer radtus of the 


graphite cell). 


Applying the condition of zero neutron current at Ps in 
Equation 11 
«a & (a oro) 
= 6 py (12) 
1 m2 
thus the moderator flux, Equation 11, becomes 
Sy 
Z.iPr) = OM, (er) + = (13) 
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where ‘ 
K, (#? .,) 


Ky (of 3) = Ky (or) + T, Po) i, (oP 3) 


Continuity of flux and current at ry gives a pair of equa- 
tions from Equations 10 and 13 which can be solved simul~ 


taneously to give 





D ww (¢#_r_) 
mh 
and G Dy # ty (of wu! co 
=~ => = (15) 
m 
where Ky (# P.) 
My Cf 3) = K, ( fr) - T, (#5) I, (oP or) 
and 


.* Z D Fo ah #0.) My #2) : Dy df. 1,4 oP.) My ( re) 


with dh (x) 


The thermal utilization, f, is the ratio of the thermal 
neutrons absorbed in the fuel to the total thermal neutrons 
ebsorbed. This ratio is aquivalent to the ratio of the ther- 
mal neutrons absorbed in the fuel to the total thermal neutron 


source assuming no net leakage from the cell. 
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th pf, &, 2a ér 
f= V So (16) 


Bi 


Combining Fquation 10 and 16 into 


i r. 
A en wu 
f= (s\ SRE Sf ty (nr) 3 ar (17) 
m m 0 


and solving the integral (1, p. 305) gives 


e= (t) 2a fy Ta Th (Mura) (18) 
"mn ™, # 


Combining Equations 1) and 18 and solving for the reciprocal 
of f gives 


1 ma + 
F = j] + — 3 FP + (& = 1) (19) 
UW 





~ oo oe 
ee - TH FT r (20) 


and 





V 
% « wm, THF, 
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F is the disadvantage factor of the urenium which is the 
ratio of the neutron flux at the rod surface to the average 


flux in the rod interior (5, pe270). 


} ir) 
Som Ay{?) “2 (22) 


U 


The average thermal neutron flux in the mederator is greater 
than the thermal neutron flux at the rod surface because the 
diffusion coefficient in the moderator is finite. The ad- 
ditional absorption is measured by the quantity (EF - 1), the 
excess absorption term. 

Murray (1, pe 100) assumed thet any peisons that ere 
tolerable do not appreciably disturb the basic neutron flux 
distribution from thet of the configuration with fuel and 
moderator only. Therefore, the aluminum cladding, coolant, 
and the process tube assembly can all be treated as poisons. 
The reciprocal thermal utilizetion with this modification can 


be written 


m1 poisons 
V, 5, 2 
1 1 
*t “oe @ (23) 
=o ta Me 
1=0 


where the average neutron flux in the added poison or absorber 


is B,- 
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Assuming that for the cladding, water annulus, and proc~ 
ess tube the ratio $,/ 4, is approximately equal to 
Z, (r)/ By Rauation 23 can be solved directly. Murray (1, 
pe 100) suggested that a slight refinement to Fouation 23 may 
be made by letting 2, be the difference between the actual 
poison cross section and the moderator cross section that 1t 
physically replaces. 

Littler and Raffle (10, pe 94) give a useful approxzima- 


tion for the excess absorption term, E - l. 





i a 3. ry 
o1w —h,f 2 . A , 
ive (~via @e toe (24) 
2 1 2 


The gonstants used in Equation 23 were listed in Tables 1 
and 2, The equivalent cell dimensions were derived from di- 
rect measurements and calculations. The thermal neutron macro- 
seopic cross sections for aluminum (cladding), water, and 
graphite were calculeted from Equation 25 using literature 


values of GO, (1), (5), and (11). 


S| 
y= Of (25) 


where 4, is 6.02 x 10°? nuclet/mole, 


and A is the mese number. 
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= for the aluminum process tube was calculated using a 
modification of Equation 25 to account for the alloying ma- 
terials in the aluminum alloy. The maxwellian distribution 
correction and a "not-1/v" correction was used to determine 
an effective absorption cross section for natural uranium. 

All the macroscopic cross section values used in this theo- 
retical analysis are shown in Table 2. 

Due to the Limitetions of simple diffusion theory, trans~ 
port theory must be used for the calculation of f° Murrey 


(1, pe 88) gives for the uranium fuel the relationship 


df of. 

—— = tanh ——4 (26) 
2 u 

x Za 


where z is the total cross section. 

Solving Equation 25 using eross section values from Ta- 
ble 2 gave a value of 0.675 en for + The inverse dif- 
fusion lengths for graphite and water were found using Fqua- 


tion 2/ and values of the diffusion lengths from the literature 


(5). 
fat= (z,/o)¥2 (27) 


The thermal utilizetionsa for the dry lattice contf'igura- 
tions were calculated by omitting the term for water in Equa- 
tion 23. The volume ratios in Iouation 19 and 23 are equiva~ 


lent to area ratios, 
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The constants listed in Tables 1 end 2 were used in 
Equations 19, 29, 23, and 2) to evaluate the comoosite ther- 
mal utilizations fer the various lattice configurations con- 
sidered. A Bessel function table (12) was used to evaluate 
the disadvantage factor, *, in Fquation 20. Equation 2) was 
used to evaluate the excess absorption term, (F - 1). 

The resulting thermal utilizations computed by Method A 
are listed in Table 3 for the various lattice confisurations. 
2. Method 3B 

A second method of calculating the thermal utilization 
of uranium is the met od of Rumsey and Volkoff (2), (3). This 
method is an extension of prior work in heterogeneous lattice 
theory (i), (6) and Includes the moderating effeet of a enter 


annulus in the lattice cell. 

















fy =f, (1 + 6) (28) 
1 
where ee FF ee + § +s, + 5, * Be 
(29) 
.¥ RN oo. V 
and 6 = tot (f= 1) + + if © t e % v 
%e “g ¥ "uw Tw “w 
" (30) 
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The “relative absorption" tern, Ry, denotes per thermal 
neutron captured in the uranium, the number of thermal neu-~- 
trons which would be captured in the jth medium if the neutron 
density in the jth medium were uniformly equal to the neutron 


density at the uranium—aluminum interface. 


os out 2. ¥ 
R, =a P= > &, 31) 
ean sh 


where F is the disadvantage factor of the uranium, Equation 
20. 

The “self blocking" term, So denotes per thermal neutron 
absorbed in the uranium, the excese number of therrml neutrons 
eaptured in the jth medium due to the exeess density of neu- 
trons in the jth medium over the neutron density at the i - jth 
interface. 

Ee ty 


Ww 


(32) 


and _* (E = 1) Liem omen en, +s, | (33) 


where (£ - 1) 18 the disadvantase factor of the moderator or 
excess absorption term given by Equetion 2h. 
The "blocking term”, Bass denotes per thermal neutron ab-« 


sorbed in uranium, the excess nunber of thermal neutrons 
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captured in the jth medium due to neutron density rise across 


the ith medium. 


N.da *V a. 
Sw Kk Pe.~ ee 2 oe 8 a2, 
Bo = OP HK & ae eee (3h) 
Wy wow By Ty “ws Ww Ww 2w iw 
nd 
a Ou o Ho. Cnr (¥ + | 
2 2 is ge al 2 can 


(35) 

In the above aix equations df, is the inverse diffusion 
length of water, and t,, is the thickness of the water annulus. 
The crosa section values are listed in Table 2, and the volume 
values per slus are listed in Table 1. The value of (7 «1) 
is given by Ecuation 2h which 4s an approximation of the com 
plex Bessel function in Fquation 21. The approximation is 
consicered to be sufficiently aceurate for the two lattice 
geometries. 

The production rate of thermal neutrons per unit volume 


per second for water and graphite are denoted Gq and q, re- 


— 


7) 
spectively. The ratio of the slowing down power of water to 


eraphite, a,/%.,» equale 20 (3, pe 22). 

Figure 1 and *ifure 2 show the neutron density éistribu~ 
tion in the dry and wet lattice configurations respectively. 
The several terms of the competitive absorption, Equation 29, 


are shown schematically in these diarrama. The placement of 
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the individual absorption terms, defined in Equations 31 
through 35, and the relative magnitude of these various terms 
ean also be seen in the dilagran.e Any absorption of neutrons 
by the air in the air annulus was neglected. The net flow of 
thermal neutrons across the boundery of a cell was asswned to 
be zero. 

B 


8 B » and 5 which are 


p’ “alw’ “alp’ “alg PE 
shown in Fisure 1 and Figure 2 were assumed asa neglisible and 


The terms Sal? 3 


were not included in Equation 29. At the uranium-aluminunm 
interface, continuity of the neutron flux density and current 
density was assumed. The neutron current density was also as-~ 
sumed es lineer through the thin mediums (al, p, or w) con- 
sidered. 

Ags an example, in deriving the Sel term, the aluminum 


surface neutron flux is 
a7 
x surface _ , | al 
Bar * fylry) * tay (*%s2) 
r 
u 
and with Equation (10) is 
D 
Surlace uU 
Bsr = (0) . sl (oe . ff. 4, ( #7) 


By definition the "self blocking" term for aluminum is 
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ave — surrecs) al 
al jb ving Vv 
u u a 


and taking from (10, p. 89) 


4 ave 2h 


and substituting sives 


x ne i 2 2 

San = Stay” Ma (37) 
Similarly 

Se = 3 Sy _. (32) 
and 


A calculation of Sai? Equation 37, showed that with #,,° = 
9.00244 em” (10, ps 93) and t., = 06102 om, that 


S = 1.27 x 107°. This value is of a magnitude which can be 


al 
considered negligible. 

The very small factor, ( f, t,)°, for the cladding and 
alw? Baip alg? ene 
_ to become negligible with respect to the other terms in 


procese tube causes the terms “1? 5 B » &B 


Figure 1 anc Pigure 2. With these eliminations, Figure 1 and 
Pigure 2 agree with Equation 29. 
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The constants listed in Tables 1 and 2 were used in 
Equations 28 throuch 35 to evaluate the wet configuration 
thermal utilizations by the Rumsey and Volkoff method. lLcua- 
tion 20 was used to evaluate * in Equation 31. The value 
(—f§ - 1) in Fevation 33 was evaluated by using Fquation 2h. 
The thermal utilizations for the dry configurations were cal- 
culated using Equatione 28 through 35 with the terms for wa- 
ter omitted. The resulting thermal utilizations, Method EB, 


are listed in Table 3 for the various eonfigurations. 


%. Resonance Escape Probability 


The fraction of fast fission neutrons that reach thermal 
energy without experiencing capture while slowins down, is 
ealled the resonance escape probability, p. The methods of 
Murray (1) and of Rumsey and Volkeff (2) were also used to 
theoretically determine the resonance escape factor. These 
methods are considered separately and are cenoted Methods A 
and & reapectively. 
1s Method A 

Murray (1, ps 92) gives the hetero, ,eneous equivalent for 


p as 





Vv # > 
psexp [- ——> re = (38) 
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where oa is the effective macroscopic cross section. 
The form of Equation 38 sugcests a resemblance to similar 
ratios appearing in the expression to evaluate f, 


A definition of the resonance utilization, f ta 


r? 


¢ = | resonance absorption in fuel 


Y pesonance absorption in fuel + seattering “ebsorption" 
oy moderator 
(39) 


The resonance absorption by the uranium is Z,, vo 2» 
from a thermal flux absorption anelogy, where 24 is an ap- 
propriately chosen average eross section of uranium over the 
entire resonance flux region. 

Assuming that the flux over the resonance region is 1/£ 


dependent, the average cross section ts 





-_ So," = C.. 
Ca” ye at. ™ i: (4.0) 
E In 


& 


where In = 5e6 for uranium metal (5, p. 273). The presence 


£ the moderator and the geometric nature of the heteroreneous 
lattice requires that an effective resonance intevral (0 W.! 


be used for the value of Oh: Equation lO is then 


~ «728 
Ou” “bet (41) 
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where (1, pe 93) 


T= S (Ty, off = 9025 + 2ke7 8/0 (42) 


re 


and S/M is the fuel surface to mass ratio in em" /g. 


oS = 
finally 24s pes (13) 


The seattering “absorption” by the moderator would be 


| Se “an = | 
mm& from a thermal flux absorption analogy, 


& gh & 
where > Fs / F z% hy di 25 (ted) 
1 SS Fab 


With Fquations 39, 43, and 4h the magnitude of the exponent 
of p in fquation 38 is 


+, Fi. one “ 1 =i 
ae = lL (45) 
Vin Aen f 2s 4 


Murray (1, p. 9h) gives an empirical expression for the ef- 


fective resonance inverse diffusion length of uranium, 
re “i 
re. * 0.0222p,, em (6) 


whore e is expressed in g/enr - 
« 
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The effective resonance inverse diffusion length for 


graphite is, 


re Mel | 
wo. pr. 2. (47) 


E 
where =e is the actual transport cross sectien of the graph- 
ite averaged over the entire resonance flux region. To 
& 
evaluate lan the relationship 


ae (1 ~ gos 9 ) (1,8) 


ie Gerived from (5, p. 98) 


A, 2 A/ = coe (9) 
where eos 6 = 2/3A (50) 


and A is the moderator mmaas rumber. 
The method of caleculatins f is now borrowed to find co 


Equation 19 modified is 


ry) © eae (52) 
rh Vea 20 | | | 


The values of wee and were determined from Equations 6 


ana 7. 2 g and 2S were calculated from Fauations 3 and kh. 
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Using the above values of ec” and to F was determined 
from Equation 20 and (FE - 1) was determined from Equation 2). 
Table 2 lista all the important lattice material constants. 
Solving Equation 51 and using the value of f, in Fque-~ 
tion kS determined the exponent of p in Equation 33. The res- 
onance escape probabilities for the various configurations, 
Method A, are listed in Table 3. The values of p for the wet 
and dry configurations were the same using Method A. 
2. Method 5 
Volkoff and Rumsey (2, p. 293) express the equation for 


the resonance escape probability as 
p=exp |- p= (52) 


where f. is the resonance utilization of uranium for reso- 
nance neutrons. Noting the form similarity of the exponent 
of Equation 52 to the therml utilization equations, the fac- 


tor, re is 


f, @ £5 (1 + 8) (53) 
with 
1 
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ehanging the cross sections on to oy and rs to ay in 
Equations 31, 33, and 35 gives a solution to Pcuation 5h. 
6 is given by Toquation 30. 

Substituting the individual terms in Equation 5h gives 


g-- lsa+y (55) 
0 
where 5 
V 
— m e, 8 ia al ecap + can 
’ Xo “ 21 “eg i Ye *w . : Ww Var 
* (BH) (2,5 Va + 2.5% 2y Vy) (56) 
and 
2, 2 1 m Ei 
Y= #°t (E) V+ v| * & p* (Bed) 
wow La. es 25 p 24 g 2 | 
(57) 
Combining Equations 52, 53, and 55 gives 
rm 1+ 6 
p = exp pore (58) 


The value of © was again obtained from Equation 2]. The 
average eross section of uranium over the resonance region, 
_ is given by Equation 3. The scattering "absorption" 
cross section of the moderator over the resonance region, a 
is given by Faquation hi. Using values frem Tables 1 and 2 
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and Equations 20, 24, 30, 43, bh, Lé, 47, 56, 57, and 58 the 
resonance escape probability was determined for Methedc B. 

These values of p are listed in Table 3 fer the lattice con- 
figurations considered. The dry valuee for p were obtained 


by eliminating the terms for water in all equations. 
P, Buekling 


The theoretical buekling was computed using Equations 1 
and 2. ‘The diffusion lemgth for each lattice configuration 


was calculated by (S, ps 280) 


rs Be aw FQ - £) (59) 
can om 
where — is the diffusion lenath in pure moderator. The val- 
we of f used in each case was the theoretically determined 
value listed in Table 3. The Fermi age, + , of the thermal 
neutrons for the lattice configurations was assumed to be the 
same as 7 for the moderator. Murray (1, pe» 123) gives the 
value of 7 as 36h om”. The value of ky, Teor eaeh configura- 
tion was determined from the values of é€ , " » £, and p liated 
in Table 3 for each con{isuration. 

Tne volume of a cubieal eritical reactor having the siven 


lattice configuration is 


vn = ) (66) 
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The sauare of the lattice diffumion length and the buckling 
for each of the lattice confirurations considered are lieted 
jn Table i. The volume and cube side longth of a cubical 
eritical reactor for each censidered lattice configuration 


are also listed in Table bh. 
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IV. EXPERIFONTAL INVESTIGATION 
A. Desesription of Equipment 


A picture of the subderitieal assembly with a 8.5 in. 
lattice spacing is shown in “igure 3. The physical dimen- 
sions of the graphite assembly were 60 in. x 60 in. x 79 in. 
The actual graphite used in the construction of the assembly 
was originally 7 in. diameter solid graphite rods. These 
rods were squared off to a 6 in. x 6 in. ecroas section (3.5 
in. radius on the corners) and stacked as shown in Figure 3. 
The top five rows of the graphite rods in the assembly were 
S in. x 6 in., and the bottom nine rows were 6 in. x 6 in. 

fhe density of the graphite was cetermined from a weight 
and dimension measurement of a mass of graphite. The graph- 
ite density was found to be 1.56 e/en’. 

The graphite assembly was covered on the top and four 
Sides with a 16 mil sheet of cadmium sandwiched between a 
375 in. thickness of plyweod and 6.125 in. of masonite. 
The cadmium sheet gave an effective "black wall" for thermal 
neutrons. Adding a 0.25 in. total spacing allowance for the 
graphite rode a horizontel section between the "black walls" 
of the assembly was 60.5 in. x 69.5 ine in area. The assem 
bly was mounted on @ concrete base, and tanks of water were 


placed under the graphite structure to insure that adequate 
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shielding was effected with the neutron source in place. 
Figure 3 shows the relative placement of these tanks. 

The lattice spacing was varied by changing the arrange- 
ment of the uranium filled procesa tubes in the assembly. 

The 6 in. lattice, Lattice I, was the conriguration with all 
the lattice holes filled with process tubes. The 3.5 in. 
lattice, Lattice II, was the confisuration as shown on Fig- 
ure 3. By removing alternate rows of the process tubes from 
the &.5 in. lattice, the 12 in. lattice, Lattice III, was 
realized. 

Five plutonium-beryllium neutron sources contained in 
individual right cylindrical containers were placed under the 
eenter of the graphite assembly. Each of these sources emit- 
ted approximately 1.63 x 10° neutrons per second; therefore 
the total source strength was approximately 6.15 x 10° new@ 
trons per second. <A cruciform geometry was chosen for the 
placement of the five sources. One source was at each of the 
cruciform ends, and one was placed at the cruciform center. 
The total source was in this way effectively contained in a 
Circle of about 3.5 in. in diameter. 

The cylindrical uranium slugs were 1.080 in. in diameter 
and 8.40 in. in length. Assuming a 0 mil cladding and «4 
206 mil slug end cap of 25 aluminum, the ectual uranium fuel 
size was a rod 1.00 in. in @iameter and 6.v0 ine in length. 


Te density of the uranium was determined from the weights 
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and dimensions of several random sample slugs. The density 
of natural uranium determined in this way was 19.9 e/ene. 

The smallest lattice spacing investigated, the 6 in. 
lattice, required 117 sections of 6061-T6 (commercial desig- 
nation 618) aluminum tubing. The tubing sections had an 
outer diameter of 14375 in., a wall thickness of 35 mils, 
and a length of 62 in. The composition of 613 aluminum is 
0.25% copper, 0.60% silicon, 1.00% magnesium, 0.25% chromi- 
um, and the remainder aluminum. 

The spacing wire composition was 25 aluminum and had a 
9.102 in. diameter. The wire was used to center the uranium 
slugs in the process tube. Seven uranium slugs were loaded 
in each tube with a helieal wrapping of wire around the slugs 
for spacing in the tube. Approximately 10 feet of wire was 
used for each seven slugs. 

Pigure shows a seale drawing of an actual unit eel 
of the subcritical assembly with a loaded process tube in 
place. Wumber seven rubber stoppers were usec at both ends 
of the process tubing to contain the water coolant. The 
stoppers were also used for the dry configurations to help 
eontain the uranium slugs during handling of the loaded proc- 
ess tubes. Allowing two process tube lengths of wire for 
each tube, the effective thickness of the process tube was 
calculated to be 0 mils. The eflecctive thickness of the 


water annulus was then found to be 3.107 in. 
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Indium foils were mounted on flat alwuinum plates which 
in turn were placed on an aluminum holder. The holder was 
designed so that the depth of penetration of each foil inte 
the slots in the assembly uld be accurately detormined, 
The indium foils were 1.5 in. x 1.0 in. x 0.093 in. and had 
en average weight of 0.5913 crams. The average thickness 
ef these foils was 6.1) mg /om™ - Pirure 5 shows the foil 
holder and foil placement during counting. The geometry of 
the foil with respect to the glass wall counting tube was 
held constant by use of the counter shelf holder shown in 
Figure 5, ) 

A model 1814 scaler manufactured by Nuclear-Chicaro and 
@ glass wall Cleger tube were used to count the beta activity 
of the irradiated indium foils. A conventional stop watch 


was used for the timing of the counting period. 
B. Procedure 


The uranium filled process tubes were loaded into the 
grephite assembly in varying arrangements for the three ceo-~ 
metric lattices investivated. The annulus between the slug 
and the tube wall was filled with water for the wet configura- 
tions. 

The neutron flux distribution was determined by the indi- 


um foil activation method. By determining the activity of a 
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foil after a long period ef irradiation and making an al- 
lLowance for decay during the procesa of counting, the satura- 


tion activity for a ,iven indium foil was obtained. 
Ago wa, e ** (61) 


where A, is the saturation activity, 
ia the activity at time t, 
A 4s the decay constant for indiun, 

and t is the time after removal from the neutron flux. 

The saturation activity is proportional to the neutron 
flux at an equilibrium state. The irradiation time of the 
indium foils for each run was in all cases at least 6.5 hours. 
This length of time corresponds to more than seven half-lives 
of indium as the half life for indium is 5h minutes. The 
indium foils used were not exactly all the same size and 
weieht, and therefore a normalising correction was made. Di- 
viding the saturation activity for a given foil by ita welizht 
resulted in a normalized saturation activity. Sinee a con- 
stant ccometry was maintained in the counter, this operation 
put all of the individual foil activities on the same basis. 
This normelized saturation activity is the activity which is 
tabulated in the Appendix, Table S through Table 16, for all 
of the foil positions in each of the lattice confirgurationa 


investigated, 
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Upen removing a foil, a three minute minimum cecay time 
was given to each foil before beginning the count. ‘This 
period of time elimineted any short lived decay components 
in the activeted indium. All of the counts taken of the 
indium foils were three minutes in duration. The decay time 
was taken as the difference between the time of the foil's 
removal from the neutron flux and tne midetime of the count 
duration. 

The foil holder slot positions on the front face of the 
subcritical assembly are shown on Figure 6. The horizontal 
x direction spacing was given letter designations A through 
J, The vertical z direction spacing was numbered vertically 
1 through 13. The actual distance in in. from the east face 
was the horizontal y direction distance designation. The 
unit of spacing was 6 in. in all cases with the sole excep- 
tion of a vertical spacing change to 5 in. at the position 
where the lattice structure changed to a 5 ing x 6 in. area. 

A glass wall Gleger tube operating at 950 volts was 
used for all or the foil activity determinations. Before 
and after each series of rungs, an operational check on the 
tube and scaler was maiie with a sample of radioactive stronti- 
um. The background activity remained essentially constant at 
an average of hl. counts per minute during the series of sount- 
ing runs. Dead time corrections were significant only on a 


few of the higher counts. 
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The surveys taken of the various lattice coutigurations 
were as follows: 
Lattices I and II, wet and dry cases. 
a. Vertical et x = 27 ins; y = 15 in. and y = 30 in. 
be Horizontal, x direction, at y = 18 in. and 
y = 30 ino; z = 30 in. 
¢. Horizontal, y direction, at x = 27 ins; 
2 = 39 in. 
Lattice III, wet and dry cases. 
@s Vertical at x * 27 tne; y = 30 ins 
Iattice without uranium. 
a. Vertical and horizontal, x direction, at 
y = 30 in. 
bd. Herisontal, y direction, at x = 27 in.; 
mz = 30 in. 
Figure 6 is a diagram showing the foil positions. 
Glasstone (5, mn. 28) showa that the flux vartetion in 
the vertical direction at distanees not too near the top of 


a rectangular subcritical assembly shape is 


giz) = ce” % (62) 


2 Lae 
y* = (=) . («) -3* (63) 
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The torn B,* fe designated aga the mterial buckling. The 
dimensions of a critical reactor can be obtained by setting 
zB“ = Be, the geometric buckling. With the material buck- 
ling equal to the geometric buckling, Be is defined as the 
buckling. The base dimensions a and b include the extrapo- 
lation length. At large distances from the assumed source 
plane the harmonic correetiona are assumed to be negligible. 

As the neutron flux is proporticnal to the indium foil 
saturation activity, a linear plot of the logarithm of 4(z) 
versus the vertical position z ylelded the quantity -7¥ 
wnich is the slops of this plot. The base dimensions were 
equal in the assembly used in this investigation. Ey edding 
an extrapolation length of 6.71 Ar. to each side of the ef- 
fective neutron "black wall", the square base dimensions be- 
came 62 in. Equation 63 can now be rewritten as 


a 


2 
B* = 2 (=) ~ ¥2 = 795 x 10°. ¥? (6);) 
where the units of Be are em". 
Rquation 6) was used to evaluate the buckling for the 
wet anc dry cases of all lattice geometries considered. The 
Slope, - § , was measured from the x ®@ 27 ins and y = 30 in. 


positions of the vertical survey plots of each configuration. 
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Ve. RESULTS 


The saturation activities for «given fcil positions of the 
conficurations considered are tebulated in Tebles 6 through 10 
of the Appendix. Plots of these values, figure 7 through Fig- 
ure 16, show the vertical and horizontal distribution of the 
neutron flux in the subcritical assembly for the Lattice con- 
firurations investigated. 

The theoretical and the experimental values of the buck- 
ling and critical reactor size for the various lattice con- 
figurations considered are listed in Table hh. Figure 17 shows 
the relationship between the theoretical and experimental 
buckling values for the wet and dry cases and the lattice spac~- 
ing. 

The y = 39 ine position of the vertical flux survey plots 
was used for the experimental determination of buckling. The 
vertical flux distribution in the s direction for the y = 18 
ine position, Figure 7 and Figure 9, was plotted only to show 
that the flux variation at this position was also exponential 
in mature. Due to the critical nature of the slope of these 
vertical flux plota with respect to the buekling value, the 
buckling from the y = 10 in. position was not caleulated. The 
buckling, Equation 6h, is a function of the square of the 
Slope, ~¥ . Therefore a very amall change in the quantity 


-¥ reflects a much larger change in the value of buckling. 
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A deviation from a true exponential nature was noted in 
the vertical flux plots at the lower and higher end positions 
of all surveys in the z direction. The lower end deviation 
was due to excess fast neutrons, 5 to 12 Mev, from the plu- 
tonium=-beryliium neutron source. All of the neutrons from 
the source in this lower portion of the assembly had not as 
yet been moderated. The neutron absorption reaction in indium 
ig for neutrons below 1.5 electron volts. The higher end 
deviation occurred approximately at the 59 in. height. This 
is the point where the lattice spacing changes for the tnaree 
different ceometric lattices considered due to the construc- 
tion of the assezbly. The slope of the flux plots increased 
near this general position in all of the vertical survey 
plots. Statistical error in counting the semll amount of 
activity at high 2 positions accounted for the relative in- 
erease in the dispersion of the counts at these positions. 
These effects can be observed on the plots shown on Figure 7 
through Figure ll. 

The vertical flux for the wet configuration wes greater 
than the vertical flux for the dry configuration up to 2 = 
36 in. for Lattice I, Figure 7, at the y = 18 in. position. 

A Similar situation existed for Lattice I, Figure &, at the 
¥ = 39 in. position with the eross over at 2 = 56 in. 

Fissure 9 shows that the wetedry crogs over point of the 
flux plot for Lattice II at y = 18 in. occurred at 2 = 30 in. 


Similarly again, the wet-dry cross over point of the vertical 
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flux distribution plot at y = 30 in. for Lattice il, Figure 
10, and Lattice III, Figure 11, occurred at 2 = he ine and 
z= 38 in, respectively. The effect of the moderation of 

the water waa creater than the effect of the water absorp- 
tion for approximately the lower half of the euberiticel as- 
sembly in all lattices investigated. The water moderation 
effect increased the thermal neutron flax by slowing up the 
fast neutrons. In the lower portions of the assembly, these 
fast neutrons are primarily from the seuree. This moderation 
effect by the water in the process tubes was more pronounced 
at the y = 30 in. position than at the y = 18 in. position 
due to the base center pesition of the source. Over the cen- 
ter, y = 30 in., the neutron flux consisted of more fast neu- 
trons at a wviven low horizental plane than other positions 
further off center. Thie excess of fast neutrons was moderat- 
@d in part by the water. At the y = 18 in. position a water 
moderation efrect was also experienced by the fast neutrons, 
but the greater graphite moderation due to the rreater dis- 
tance from the source caused the number of fast neutrons to 
be smaller than at the center positions. The number of fast 
neutrons due to fission was also larger at the center poai- 
tions than at the outer positions of the assembly. This is 
Gue to the sine distribution of the thermal neutron fiux in a 
horizontal plane in the asnembly. Thernal neutrons caused 
the fissions in the natural uranium which produced fast nen-~ 


trons and other fission products. The position of the maximum 
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thermal noutron flux was therefore also the position of the 
maximum fast neutron flux due to fissicn. 

The vertical flux survey for the case of the graphite 
lattice without uranium is shown on Figure 12. The dry con- 
figuration flux values for Lattices I, TI, amd III were cross 
plotted to indicate the cross over point on the flux plot 
without uranium. The flux plot of Tattice ITI crossed the no 
uranium lattice plot at a slightly higher point than the flux 
plot of Lattice II did. The flux of Lattice I head the lowest 
cross over point of the three lattices. The increases of neu- 
tron flux due te subcritical multiplication in fattice I ap~« 
peared to be leas than the absorption effect of the uranium. 
This effect was not nearly so pronounced in Lattices II and 
IfI. The neutron flux plot for the no uraniwn lattice ap- 
pearec to be truly exponential in character for all but the 
low and high extreme positions in the assembly. 

The slopes of the vertieal flux plots for the x = 27 in. 
and y = 30 in. position were determined, and it was found 
that the marnitude of these slopes varied as follows: the 
flux plot in Lattice II had the smallest slope and the flux 
plot in Lattice III had the greatest slope. The slope of the 
flux plot in Lattice I was slightly greater than the Lattice 
iiI slope. The greater the slope of the flux plot, the smal- 
ler the value of the buckling becomes. This slope varietion 
wes trie for both the wet and dry conivigurations and wes di- 


rectly reflected in the buekling values calculated. 
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There was little difference in the effect of suberitical 
multiplication betwoen the dry cases of [attice IZ and ITI. 
From a vertical height of kS in. and upwards, the flux in 
éry Lattice II wae slightly higher than in dry Lattice III. 
The absorption difference in these two lettices was apparent- 
ly small. 

The horizontal surveys for symmetry in the x direction 
for z = 30 in. are plotted in Figure 13 through Figure 15. 

An extrapolation length was added to each end beyond the ef- 
fective neutron “black wall“. All of the symmetry plots were 
plotted as sine curves. 

The horizontal x direction Lattice I flux, Figure 13, for 
the wet case was ecreater than the flux for the dry case at 
both y = 18 in. and y = 30 in. positions. The horizontal x 
direction lattice Il flux, Pigure 13, showed that there was 
very little difference between the wet and the dry cases rer 
either the y = 18 in. or y = 30 ins pesitions, The y = 30 in. 
position flux for both lattices was creater in magnitude than 
the y = 16 in. position flux due to the sine distribution in 
the y direction. 

The effect of the water moderation appeared bo be greater 
than the water absorption effect in Lattice I at both the 
y = 18 in. and y = 30 in. positions for the x direetion survey. 
in the x direction survey for Lattice II, these effects were 


epproximately equal for both y positions. 
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Lattiee II flux in the x direetion at 2 = 30 in. had in 
all cases almost en ideal sine distribution. The flux in 
wet Lattice I had some random scatter over the center sestion 
f the horizontal survey. This variation pernaps was cue to 
water moderation having as preater effeat than the water ab- 
serption over the near center section. This net effect was 
more pronounced directly over the sources. 

Pigure 1h shows a horizontal x direetion cemparison at 
2 = 30 ine of the wet, dry, and without uranium confirurations. 
teattice I values of flux for both the wet and dry cases were 
below the flux values in the no uraniun lattice, The flux in 
wet Gattice I was again higher than the flux in dry Lattice 
is In fattice I, the absorption effect at the « = 30 in. 
position was obviously greater than the suberitical multi-~- 
plication efrect. Lattice II values of flux for the wet and 
ary cases are Just alightly below the no uranium case. In 
this lattice the suberitical multiplication effect and the ab- 
sorption effects were approximately equal. 

A lattice comparison for the horizontal x direction at 
2 = 30 in. is shown on Figure 15. The Lattioe II flux was 
greater than the Lattice I flux in both the wet and dry con- 
figurations. 

Symmetry (gine distribution) was shown in the y direc- 
tion on Figure 16 for Imttices I and II at the x = 27 in. and 
® # 30 ins position. The plet shows that there was little 


difference between the wet and dry cases of either Lattice I 
















OW ok 0 ee rer © rd i endl OC oandinne 
8 SF SER cet) Jems ce ately 
A122 SO =H EA TS. CBs hed § ohhbde) Oe 
2k cee, teeter ae (pers 
ee ee ee 
oe Mie cme Oh cele wine eee et oe 
: a awe a hee ths oem 
1) CPA ke tlhe) « dake Ge 6 cmt 0F omrel's 
ite) 2 eee Gs ae fF ok a Te 
me ences GS tee Ce et me er eee fe epeter 4 
oe ee ae, | 
a Ted ele th nt ieee sees com § onion 
(Oh Oh @ de yee eerkeemede wee Ue weld ond oe 
he hag om a tow Chere somes cp 
= ee at teeree O sateen kee 
ee. eee ee 
0 BO ce me te em ee hk te oe (ebm ole mde der @ 
Os 0 el ee 
oa) —_hm a! et miiemm, amritel + 
wend ees ae OS ree oe oe ol 8 
ee yer? Ge tee oe! eee «) ant” 6 eee et eek 













~—_ii 








é 


eT Oe ok cee ee (eee diet onze) autem 
ee Oe ee 
DOM ot 1 EN Ae a MMe ot Oe 
Tr. caer © lee? er hee OP Ge Ag) qummretoNe 





UT 


or II. In tattice I there was an apparent digression from 
a sine distribution at the center positions of the wet con- 
figuration. Also the sine curve plotted for symmetry was 
slightly to the left of the peints pletted on Figure 16. 
This was probably due to the placement of the foil holders 
or gome spreading of the graphite blocks from a Lattice I 
loading. The flux plots of the wet and dry configurations 
in the y direction for Lattice IIT are approximately the same 
as the flux plot of the no urenium lattice. There is asain 
some dispersion at the center position in this comparison. 

The results of the experimental and theoretical bucklings 
determined for all configuretions considered are shown on Fig- 
ure 17, The two theoretical mothods used to determine the 
buckling for the dry confisuretions showed reasonable arree= 
ment, but the wet case theory was in disagreement. The re- 
sults from Method A were hifh in comparison with the results 
[rom Method B for wet Lattices I, II and III. 

The experimental buckling values for the wet configure- 
tions were lower than the experimental buckling values for 
the dry configurations in all three lattice spacings investi- 
feted. The theoretical values of buckling determined for the 
dry configurations, Methods A and B, were a little higher 
than the experimental values of buckling for the dry confirura- 
tions. The theoretiesl values of buckling obtained from 
Method 4 for the wet configurations were high in eomparison 


with the experimental buckling values. The theoretical Method 





















(A) (GIG IS HY HA Wee me! OWE OD. 74 
ee ee TT ‘a 
a TT mh ee come ent wate 
(O! ate «thee e nko oom bem bbe onde of © 
—t tt a 1. Aes a oo bok chee 
+ ie FF ne men. Foe amnowes 
ee en eee le 
= oO (eter = «te |” en rem aL atlenenet ¢. 
lee 60 ee les eee ee ee tml wel 
_ AR ht AON Nhe. ae th Nn enermaye 
ett! gmat = sere ems «ee % 6Meers a 
i ele eee ; 
oF ee tee ete leetdeeee are e@® LE 
ee eee ee ide me Oe oe ae ett : 
St lel de nd ee bee ee ane 
is FP et eee © 6 tt om 6 ete ah 
« - eueme G 4 aieue ae wh © meme. 
— es ww =< = ei tte (Aterpemes aa 
ee) et dees tte 6 Ge care comm) eoom 
ttt eh debe! oe Vw oe eedeemenh ome ooh 
7) Ce ime © eee eee Ae 
wee cee) ogee .) © ' (eee 
7 ee ete || ee lemme ane 
ee ee oe 
eee lL eb come cet) om lee fee tr ~~) « 
Ci le nn he) 
























4.8 


B gave buckling values for the wet confiscurations which were 
in excellent agreement with the experimental buckling values 
for the wet lattices. 

For the natural wraniuma graphite assembly a maxinun did 
oceur for the wet and dry cases of the buckling versus lattice 
spacing curve, Pigure 17. This maximum cccurred in both the 
experimental and theoretical analyses. The C/U aton ratios 
were 73 for Lattice I, 16 for Lattice II, and 292 fer Lattice 
III. Davenport (8, p. 315} coneluded thet the maximum bucke 
ling for e given natural uranium fuel assembly should occur 
at a value of ¢/U between 50 end 100. The maximun buckling: 
in this investigation of three lattice spacings occurred at 
about a C/U value of 16 for both the wet and dry configura- 


tions. 





9 


VI. DISCUSSION OF RESULTS 


The theoretical methods employed had the following prob- 
able errors. Accurate detailed neutron scattering and capture 
eross section data are not available in the literature for all 
materials. Corrections for varying densities and purities of 
the materials were not computed. A calculated disadvantage 
factor for uranium was used instead of a measured value. Val-~ 
ues from the literature were used for the Permi age and the 
various lattice diffusion lengths where experimental values 
again would have given more accurate results for the calcu- 
lated buckling. There were cracks between the graphite blocks 
with a proeess tube uranium load in the assembly. These air 
spaces and the air effect around the process tube itself were 
assumed to be negligible. The effect of temperature and hu- 
midity on the results was also assumed to be negligible. Fi- 
nally, one group neutron theory was used throughout the analy- 
sis. 

The following comments can be made concerning the experi- 
mental methods employed. The harmonic corrections to the flux 
distributions obtained were assumed to be negligible. The 
positioning of the foils and the normalizing of the foil ac- 
tivity were subject to human error. Finally, and most impor- 
tant, the slope of the in @(z) versus z direction curve must 
be determined with great accuracy. This slope value is very 


critical in determining the actual buckling. Repeat runs for 
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consistancy et a given position should te taken from which 
an average value of activities at a position could be used 
for plotting ard the POR Any SAGA Sie, Varyiny; the foil 
weight at a given position during the runs would tend to can- 
cel out any normalizing error. 

The feasibility of an actual reactor with the process 
tubing, natural uranium, slug geometry, and water coolant used 
in this investigation was determined. Assuming 1250 channels, 
a water coolant, and 0.105 in.* flow area per channel the mass 
flow is 1.68 x 10° 1b /hr. Aasswning a temperature difference 
of 70°F and using t, = 1,06 for water at hOO°F, the heat rate 
transfer is 366 megawatts. This heat rate transfer is large 
enough for a potential power reactor. The above assumptions 
are conservative in comparison with the experimental critical 
reactor volume determined for the wet 8.5 in. lattice config- 
uration. 

The process tubing did spread the graphite blocks a small 
amount due to the fact that some of the holes were smaller 
than the normal average designed size. A power drill of an 
appropriate size could onlarge these holos so that the process 
tubing would slide in and out of the assembly with ease and 
ne spread would occur. It was found that a lattice change was 
more easily accomplished using the process tube uranium con- 
tainment than with just the loose slugs. One end of the proc- 
ess tube was beveled to facilitate loading of the wire wrapped 


Slugs. It is sugcested that a number seven rubber stopper be 
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well seated in the beveled end and a larger stopper be used 
at the other end. Coolant leakage would in this way be elimt- 
nated. 

The E=3 pesition foil holder hole was a through hole in 
the y direction. This through hole would be of more use if 
it were moved up to the £-8 or F-9 positions. The H=3 posi- 
tion was too close to the source for 4a food symmetry flux 
survey. 

The 6.5 in. lattice is the optimum lattice spacing which 
is possible with the subcritical assembly used. Suggested 
ideas for future study are the investigation of other coolants 
and a more detailed enalysis of the §&.5 in. lattice for both 
the wet and dry configurations. Also, the effect of the re- 
flection by the water contained below the source may be in~ 
vestigated to cive a maximum source flux with the evailable 
sources in place. With the use of an additional inner process 
tube the effect of a variable wator annulus could also be in- 


vesticated. 
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Vits CORCLUSIONS 


The moderation effect of the water was creater than the 
absorption effect of the water Cor approximately the lower 
half of the subcritical assembly for the three lattice spac- 
ings investigated. The thermal neutron flux in the lower 
portion of the assembly was grecter in the wet con‘igurations 
than in the dry configurations, and the reverse situation 
existed in the upper portion of the assembly. The 6 in. wet 
lattice hac the loweat overall flux because the absorption 
effect of the greater amount of water and uranium was larger 
than the subcritical multiplication effect in thia lettice. 
The 12 in. lattice was "over moderated" with an apparent small 
suberitical multiplication effect. The optimum lattice in- 
vestigated was the 8.5 in. lattices. In this lattice the water 
had very little effect on the Tlux distribution due to ths 
balance of the water moderation and absorption effects. 

The symmetry surveys showed that in the x and y direc- 
tions that the flux did follew a sine distribution for both 
the wet and dry configurations. 

The two theoretical methods used to determine the buck- 
ling for the dry configurations fave results with reasonable 
agreement. These theoretical bucklings for the dry cases were 
high in comparison with the experimental bucklings determined 
for the dry cases except for the 12 in. lattice. The two 


theoretical methods used to determine the buckling for the 
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wet configurations gave results which were in disagreement. 
The method of Murray gave values o: the buckling for the wet 
configuration of the 6 in. and 8.5 in.lattices which were 
very high in comparison with the experimental buckling values. 
The method of Rumsey and Volkoff gave values cf the buck-= 
ling for the wet confisuration of the 6 in. and 8.5 in. lat- 
tices which were in excellent agreement with the experimental 
buckling values determined. The experimental buckling value 
determined for the 12 in. wet lattice was between the two 
theoretical buckling values determined for this lattice. The 
experimental buckling values for the dry cases were higher 
than the values for the wet cases in 211 lettice spacings 
investigated. The only experimental positive value of buck- 
ling for the wet configurations was the buckling value for 
the 8.5 in. lattice. 

A maximum in the buckling versus lattice spacings curve 
did occur for both the wet and dry conficurations. This 
maximum in buckling was found in both the theoretical and the 
experimental values determined. These buckling curves for 
the wet anc dry configurations did not cross; hence a "fail- 
safe" behavior upon the loss of water coolant could not be 


expected from these results. 
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Lattice I (6 ins), II (8.5 ine), and Iii (12 in.) 
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Table 1 


Dimensions for a unit cell 








Ps uranium rod radius 

toa » thiekness of aluminum slug can 

t 5 effective thickness of water annulus 

t,, > effective thickness of process tube 

veir , effective thickness of air annulus 

ri equivalent inner radius of graphite 
Foquivalent cross section areas for one cell 

Uranium 

Slug cen 


Water coolant 


Process tube 


A 


ir 


Graphite, imttice TI 


Graphite, Lattice II 


Graphite, Lattice IiI 


Volume per slugs 


<g 


u 9 


uranium 

slug can and cap 
wator 

procesa tube 
graphite, Lattice I 
erephite, imttice Ii 
graphite, Lattice III 





ae oi Ae 


1.270 om 
0.192 om 
e273 em 
O.102 cm 
06455 om 
20620 cm 
5.2% em 
0.8 em@ 
2.61 om 
1606 om® 
5.65 om* 
217.0 om 
434s om™ 
668.0 em 
103.9 om 
23.0 ew 
55.7 one 
22.6 em 
hé2s5 en 
9250 ene 


18,500 em? 
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Table 1 (Continued) 





Lattice I, rj, equivalent outer radius of graphite 8.60 om 


Lattice II, r., equivalent outer radius of graphite 11.9) em 


2 


Lattice Iil, Pos equivalent outer radius of 
grephite 16.76 em 
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Table 2 


Lattice cell material constants 


Symbol Value Sour ee 

ri 0.01323 em™+ (5) and Equation 25 

a 0.01556 om™? (5), (11) and Equation 25 
r 0.32). en” (1) and Equation 25 
- 0.399 em” > (1) and Equation 25 

Zn 0.0967 om Equation 13 

xe 0.00036 em + (5) 

p> 0.375 em”! (1) and Equation 25 

>; 4,01057 em” * (1) and Feuation hh 

>t 0.354 em”? (1) and Zequstion h9 

z=, 0.017 om (S) 

#, 0.675 en”? Equation 26 


d 9.01992 om (5) and Pouetion 27 
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fable 2 (Cont inued) 





Symbol Value Souroe 
of. 063472 em™* (5) ané@ Pquation 27 
Ly 0.8218 en” * Equation 6 

df - 961059 em”* Equation 47 
cos 9 0.0555 Fquation 50 
a,//4, 20 (3) 

F. 0.158 (5) 
Fermi age, v. 36h en (1) 

ue 2520 em@ (1) 

(. 1.256 fom Measured value 
Cs 19.0 g/one Measured value 


s/i «0628 em®/em Calculated value, Equation 2 





= «a 
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*y valuated using €é= 16927 and 7 = 1.308, 


6926 
0.823 
1.923 


0.901 
0.823 
0.996 


929 
0.836 
1.03 


0.682 
0.832 
0.985 
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Table 3 


Method A 


Method B 


Theoretical nuclear lattice conatents 


A pa TURE ATES A a a 


sms 


0.870 
06909 
1.062 


0.67 
6.909 
1.035 


If 


6. &63 
0.915 
1.061 


06818 
0.911 
1.008 





5.767 
6.955 
0.984 


9.751 
00955 
0.963 


IIT 


0.759 
0.958 
0.976 


6.710 
6919 
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Table 
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Theoretical Method A 











Lattice I il iil 
DRY Value 
1? om® 186 328 587 
p “en hl x 107° 86 x 107° «17 x 107° 
Ven? 611 x 10° 20: x 10° - 
Cube side, ft. 27te7 L$ ~ 
wit? Value 
Len* 250 386 627 
BR *em™* ~6.5 x 1076 LO x 107° =37 x 107° 
Ven? - 620 x 10° - 
tn de sll ERE ee i+ nt &- 
Theoretical Method 8 
Lattice L It iit 
DRY Value 
tL em® 201 3h5 607 
Bon 59 x 107° 63 x 107° -25 x 107° 
Ven 357 x 10° 212 x 10° = 
Cube side, ft. 23.62 19.6 ~ 
WELT Value 
ten 298 LSg 730 
fon” “23 % 1076 9.7 x 1076 ~O5 x 1976 
V nen - 530 x 10° ~ 


Cube side, ft. - 57k - 
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Table |. (Continued) 








Experimental 
Lattice I It Tit 
DRY Value 
Bom“ 35 x 1976 70 x 1076 “5 x 1076 
V,00 778 x 10° 275 x 10° ~ 
Cube side, ft. 30.1 21.3 ? 
WET Value 
Ben “23 x 1076 l2 x 107° #8 x 107° 
V,on° “ 3680 =x 10° - 


Cube side, ft. - 51.5 - 
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Figure 3. Suberitical assembly 
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Figure 6. Front view of assembly showing foil positions. 
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Figure 7. Vertical flux survey, Lattice I, 


x = 27 ine, y = 18 in. 
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Figure 8. Vertical flux survey, Lattice I, 
x = 27 ine, y = 30 in. 
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Figure 9. Vertical flux survey, Lattice II, 
x= 27 ine, ¥ = 18 in. 
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Figure 10. Vertical flux survey, Lattice II, 


x = 27 ine, ¥ = 36 in. 
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Figure 11. Vertical flux survey, Lattice III, 
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Figure l2. Vertical flux survey comparison, 


x = 27 ine, y = 30 in. 
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Figure 13. Horizontal flux survey, x direction, Lattice I and 
Il, zz >= 39 ine, J = 18 in. and 7 = 30 ine 
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Pigure lhe Horizontal survey, x direction, configuration 


comparison, y = 2 = 30 in. 
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Figure 16. Horizontal survey, y direction, configuration 


comparison, x = 27 ina, 2 = 30 in. 
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G/é c/i 
Position DRY WET Position DRY WET 
E1-18 3903 350 F1=30 6508 7189 
2-18 30k0 3395 E2=-30 4560 S110 
E3-18 2192 24609 E330 3176 3493 
Fh=18 1580 1661 a3 0 2085 2256 
E5~18 1029 1095 E5=30 1296 1482 
26-18 716 707 E6=30 888 870 
E7-18 65 1.78 E7=30 552 G20 
B6-18 335 313 BE=30 376 378 
£9-18 238 239 B9-30 268 253 
FL0-18 135 115 F10=30 163 158 
B11-18 107 91 F11-39 129 Lhé 
FL2-18 59 62 E12=30 69 7 
E13-18 28 32 13-30 bd 56 
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Horizontal flux survey Lattice I 
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Position DRY Wit Position DRY wet 
A5=-15 200 226 A5=30 256 321 
B5-16 89 469 B5-30 6605 650 
C5-18 678 765 C5=30 956 990 
nS-18 &93 950 D5=30 1223 1237 
E518 1029 1095 E530 1290 182 
P5-18 1030 121k5 F5=39 133% 1453 
G5-18 952 1929 G5=<30 1163 1270 
HS=-18 795 728 H5=30 939 99k, 
I5~18 505 Soh 15=30 623 671 
I5=18 21h 223 I5<30 272 282 
B5<9 115 73 £5=36 1306 Lhé5 
¥S=6 h12 yds? ES-h2 1131 1131 
ES-12 759 709 £5“eh6 875 91h 
FS-18 1029 1095 Eo=-Sh 525 583 
EB5=—2h, Lah? 1320 ES5=60 180 150 
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Table 7 


Vertical flux survey fattice II 





Saturation Activity Saturation Activity 
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eat rae warty 


c/s aft 
Position DRY WET Positio:: DRY WET 
£1-18 4S90 h9OO E1=39 700 7815 
E2=-18 3665 3802 R2-30 5505 5690 
E3-18 2720 2790 3-30 3768 3975 
E-16 1899 1930 Theo 2689 2607 
E5<18 12h 1235 L5-30 1583 1571 
E6<18 BOL 78h 6-30 1009 10hg 
E7+18 536 499 e730 707 651 
E8-18 376 332 £8<3% hho bS3 
E9-18 255 208 E9=-30 282 278 
£10=16 135 147 B10=39 175 186 
E1118 111 gé E11+36 135 118 
E12-18 79 62 B12-30 100 96 
E13-18 35 38 £13430 50 Lé 
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Horizontal flux survey Lattice II 


Saturation Activity Saturation Activity 
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Position DRY . Wit Position DRY Wet 
A5=-18 262 252 AS~30 3h0 30 
D518 632 607 BS=39 791 728 
c5=18 898 926 C530 1170 1169 
DS<18 1265 1055 D5@30 1456 1435 
E5=-18 W2k1 1235 ES-36 1583 1571 
r5<18 1239 1227 P5=39 1578 1550 
G5-18 1111 1156 G5—=3) 11,87 Lhé0 
ES=18 891 886 5-39 1129 1139 
516 602 595 15-30 762 759 
J5=18 266 221 T5=36 308 291 
E5<0 117 11:8 E5=36 16h1 1572 
ES-6 511 Sul B5<12 137 1360 
wG-12 G87 G7 BS—46 965 1002 
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B5-2h 1550 E560 195 160 
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Table 9 


Vertical flux survey Lattice III 


SP BORNE LP, 





Saturation Activity 


c/M 
Position DRY WET 
£1-30 7500 7985 
E2-30 59h5 6150 
£3<30 3960 3999 
Eh-30 2813 272 
E5=30 1660 1658 
B6=30 1152 1093 
E7=#30 669 637 
B8<36 458 h2s 
9-30 265 23h, 
F10-30 18h. 153 
#11-30 148 118 
E12=30 19 51 
E1330 9 56 
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Table L192 


Lettice flux survey without uranium 


Position Saturation Activity 
c/™ 
B1-30 8365 
2-30 6290 
B3-39 ulo8s 
Ely-3 0 2795 
F5=30 1655 
F6-39 965 
E7=30 587 
£5=35 361 
E9=30 176 
B10=30 128 
11-30 BL 
Bl2=30 Sh 
13-30 22 
A5=30 320 
B5—30 736 
05=30 1161 
DS5<3C 150 
¥5=39 1669 


Position 


"= 


£5=-9 

ES<6 

BS-12 
E5<18 
ESm2y 
£5=36 
BE5-L2 
ES<h8 
E5-5h 
E5=60 


G5=36 
H5=30 
I5=30 
J5=30 


Saturation Activity 


C/¥ 
172 
566 
992 
1339 
1590 
1639 
1320 
990 
615 
183 


1h86 
113 
735 
318 
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